Bound states of two bosons in an optical lattice near an association resonance 
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We model two bosons in an optical lattice near a Feshbach or photoassociation resonance, focusing 
on the Bose-Hubbard model in one dimension. Whereas the usual atoms-only theory with a tunable 
scattering length yields one bound state for a molecular dimer for either an attractive or repulsive 
atom-atom interaction, for a sufficiently small direct background interaction between the atoms a 
two-channel atom-molecule theory may give two bound states that represent attractively and repul- 
sively bound dimers occurring simultaneously. Such unusual molecular physics may be observable 
for an atom-molecule coupling strength comparable to the width of the dissociation continuum of 
the lattice dimer, which is the case, for instance, with narrow Feshbach resonances in Na, 87 Rb, and 
133 Cs or low-intensity photoassociation in 174 Yb. 
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An optical lattice [1] modifies the motion of atoms pro- 
foundly compared to free space. Moreover, the site-to- 
site hopping of the atoms can be controlled by varying 
the intensity of the lattice laser, which enables phenom- 
ena ranging from Mott insulators [2] to atomtronics [3]. 
The interactions between the atoms can also be tuned 
by utilizing Feshbach [4, 5] and photoassociation [6, 7] 
resonances. In effect, both the mass of the atoms and 
the atom-atom interactions can be controlled experimen- 
tally, possibly leading to custom-tailored molecules and 
unprecedented experimental control of molecular physics. 
As the prime example to date, dimers bound by repulsive 
atom-atom interactions have been demonstrated experi- 
mentally [8] and confirmed theoretically [8-12]. 

In the conventional single-channel description of the 
lattice dimers [8, 10-16] the atoms have an interaction be- 
tween them characterized by a scattering length that di- 
verges at the resonance. The more nuanced two-channel 
theory [9, 17, 18] asserts that there are also molecules 
present as an independent degree of freedom. In this 
view atom pairs may be converted into molecules, and a 
resonance occurs when the atom pairs and the molecules 
have the same energy. 

In the present communication we formulate and solve 
the time- independent Schrodinger equation for a lattice 
dimer within the two-channel model. We find a pecu- 
liar qualitative change in the molecular physics com- 
pared to the single-channel model: while the latter al- 
ways presents one bound state for the dimer, in the case 
of sufficiently weak direct interactions between the atoms 
the two-channel model may exhibit two bound states, 
one below and one above the band of dissociated dimer 
states. In particular, while the single-channel description 
provides for a bound state below the continuum band 
for attractive atom-atom interactions or a bound state 
above the dissociation continuum for repulsive interac- 
tions, the two bound states for the two-channel dimer 



may be viewed under certain conditions as the analogs 
of attractively and repulsively bound pairs occurring si- 
multaneously for the same system parameters. 

We begin with the Bose-Hubbard Hamiltonian for a 
one-dimensional lattice in the single-channel picture [11], 
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The index k runs over the lattice sites, L of them. We 
use periodic boundary conditions, so that the site k = L 
is the same as the site k = 0. The parameter J charac- 
terizes tunneling from site to site, U quantifies the atom- 
atom interactions, and a k is the annihilation operator for 
bosonic atoms at the site k. 

The corresponding two-channel model describes an as- 
sociation resonance where a bound molecular state and 
a state of two asymptotically free atoms become degen- 
erate. The Hamiltonian is 
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where b k is the annihilation operator for a molecule at 
the site k, and the detuning 8 quantifies the difference 
in energy between a molecule and an atom pair; 6 = 
denotes the resonance. The entities annihilated by a k and 
bk are bare atoms and molecules that would exist without 
the atom-molecule conversion ex £. Diagonalization of 
the Hamiltonian (2) gives a description of the physically 
observable dressed molecules. 

In our simplest possible model we do not include site- 
to-site tunneling of the molecules. This is reasonable 
since the molecules are twice as heavy as the atoms and 
the tunneling amplitude J is exponentially small in mass. 
There could also be atom-atom interactions as in the 
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single-channel model (1), as well as interactions between 
molecules and between atoms and molecules, but these 
are all ignored in the Hamiltonian (2). In Eq. (2), and in 
most of this Rapid Communication, we assume that the 
resonance at S ~ dominates the physics to the extent 
that direct two-particle interactions can be ignored. 

Although the single-channel model speaks of only 
atoms that interact among themselves and the two- 
channel model shows atom-molecule conversion with 
seemingly no atom-atom interactions at all, there is 
a close connection between the two descriptions. To 
demonstrate this, we write the Heisenberg equations of 
motion for and bk from the Hamiltonian (2), and for 
an asymptotically large detuning S eliminate adiabati- 
cally the molecular operators bk- With the identification 
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the result is the same atomic dynamics as per the single- 
channel Hamiltonian (1). One may think of the single- 
channel model as the limit of the two-channel model far 
away from the resonance. 

The analysis of the two-channel model proceeds much 
along the lines of the single-channel case [11]; the mul- 
titude of mathematical complications is essentially the 
same, and we only outline the main points. First, us- 
ing the discrete Fourier transformation we take both the 
atomic and the molecular operators into the lattice mo- 
mentum representation, — > c q and bk — > d q . The 
quasimomentum q runs over a suitable set of values, e.g., 
q = 2im/L for integers n such that q belongs to the first 
Brillouin zone of the lattice. The Hamiltonian becomes 
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For diatomic molecules, the total atom number 
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is a conserved quantity for the two-channel Hamiltonian. 

Given two lattice momenta q\ and q2, we next define 
the analogs of the center-of-mass (cm.) and relative mo- 
menta P = qi + q2 and q = (qi — 52) /2, and write an 
ansatz for the state with the total number of atoms N = 2 
as 
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Here |0) is the particle vacuum, and without restricting 
the generality we set A(q) = A(—q). The Hamiltonian 
maps the state (6) to a state of the same form, with 
the same P, which is a manifestation of conservation of 



the cm. momentum. The problems with the notation 
± q [11] are irrelevant in our final limit when the 
momentum becomes a continuous variable. Defining a 
characteristic frequency f2p and scaling all dimensional 
quantities to it, 

A- (7) 

we find the time independent Schrodinger equation for 
the energy E <-> uj in the form 
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E(3 



2L 



0. 



(A-w)0-^Y, A M = O - 
* 1 



(8) 
(9) 



From the Schrodinger equations (8) and (9) one may 
deduce the condition for the (scaled) eigenstate energy 
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For comparison, the single- channel result is 
f(uj,L) = l/lC, 



(10) 
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with K, — U/Vtp. Once the eigenenergy w has been 
solved from Eq. (10), the unit-normalized solution to the 
Schrodinger equation is 
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We plot the function f(uj, L) for L = 16 as a func- 
tion of the variable uj in Fig. 1. The right-hand sides 
of Eqs. (10) and (11) for the parameters K, = A = 1/2 
and S = l/\/2 are also represented as straight lines in 
Fig. 1. The solutions occur where the straight lines and 
the graph of f(uj,L) intersect. Equations (10) and (11) 
both have one solution uj between each value of — cos q 
for the successive discrete values of the quasimomentum 
q. This is the finite-lattice analog of the dissociation 
continuum of the dimer. The single-channel eigenvalue 
equation (11) also has one solution uj outside of the range 
(—1,1). For a positive atom- atom interaction coefficient 
U this is the celebrated repulsively bound dimer [8]. On 
the other hand, the two-channel model as in Eq. (10) al- 
ways has two bound states; one above and one below the 
dissociation continuum. 

With the addition of the explicit molecular degree of 
freedom, for each fixed cm. momentum P the dimer sys- 
tem has one more degree of freedom in the two-channel 
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FIG. 1. (Color online) The function f{u,L) [Eq. (10)] for 
L — 16. The horizontal and positive-slope lines represent 
the right-hand sides of Eqs. (11) and (10), respectively, for 
K. = A = | and 3 = l/v2. The dashed vertical lines are the 
asymptotes of f(ui,L) at the values of uj = — cosg such that 
f(u),L) = ±oo. 



case than in the single-channel case, which accounts for 
the existence of the extra state. However, the nature of 
the newly emerged bound state calls for an investigation. 

We discuss the limit when the lattice is infinitely long, 
L — > co, whereupon the quasimomentum becomes a con- 
tinuous variable in the interval [— 7r, tt). The main tech- 
nical rule is that any sum over the lattice momenta is 
replaced with an integral, ^2 — > ^ dq . Here we 
only look into the two bound states with |w| > 1. The 
eigenvalue equation and the amplitude of the molecules 
in the bound state are 
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References [9] and [18] based on scattering theory also 
report multiple bound states, but the explicit results are 
different from ours as these papers incorporate a signif- 
icant "background" atom-atom interaction cx U as in 
Eq. (1). We may add the background interaction into the 
present two-channel model, and represent it in terms of 
the same parameter JC = U/flp that was used in Rcf. [11]. 
For instance, the counterpart of Eq. (14) then becomes 



sgn (lo) 

K[u>- A) + S 2 _ V^ 2- ! 
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The advantage of the present approach is that, once we 
have both the bound states and the dissociated states, 
we can straightforwardly analyze [11] quantities such as 
the dissociation rate of the bound state when the lattice 
parameters arc modulated, as measured in Rcf. [8]. How- 
ever, the continuum states present similar mathematical 



issues as in Ref [11], and we defer the discussion of the 
dissociation rates for later. 

The variety of special cases that may arise with JC ^ 
complicates the discussion, so we continue without the 
background interactions and set JC — 0. Given the actual 
parameter values, the number of real roots to Eq. (16) 
with | a; | > 1 may be employed as a test of the qualitative 
validity of the present arguments. For instance, on exact 
resonance A = there arc two such roots if S 2 > \K\. 

We approximate the bound states as a function of the 
scaled detuning A. First, we note that if the pair (uj, A) 
satisfies Eq. (14), then so docs (— lo, — A). In view of this 
invariance, it is convenient to assume first that A > 0, 
and then use the symmetry to deal with the case A < 0. 
In the limit A 4 oo the angled straight line in Fig. 1 
slides down, and it is evident that the two bound-state 
energies become lo ~ A and lo ~ — 1. Let us first study 



the state with 



A. Since an expansion in A 



is 



expected in the limit A — > co, we attempt to find the 
energy by substituting an expansion of the form u>a = 
A + a + a_iA _1 + . . . into Eq. (14). MATHEMATICA 
makes this sort of work extremely easy. The result is 
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An expansion of the other bound state energy lot — — 1 
that approaches the continuum threshold in the limit 
A — > co is found similarly. Analytical expansions are 
also available in the neighborhood of the association res- 
onance A ~ 0. In this case the energies ui± are found in 

the form oj± = ± sgn(w) &o + b±A + 

The extra bound state can be characterized starting 
from the asymptotic expansions of the bound-state en- 
ergies. For large detuning the expansions (17) and its 
counterpart for the energy lot give the molecular frac- 
tions / = |^| 2 
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whereas for A 
states are 



the molecular fractions for both bound 



f± 



2E 



1 + 4S 4 + VI + 4S 4 



(19) 



According to Eq. (3), the single-channel model is re- 
captured from the two-channel model in the limit when 
the parameters A and 3 both tend to infinity in such a 
way that fC = — 2 2 /A remains constant. K. is then noth- 
ing but the two-channel counterpart of the dimensionless 
atom-atom interaction constant that we also denoted by 
K, in our single-channel theory [11]. In fact, in this partic- 
ular limit the energy of the bound state lot converges to 
the single-channel result lo = sgn(/C) yl + JC 2 . Moreover, 
the molecular fraction fx of this bound state vanishes. 
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The |A| — > oo bound state that we have denoted by the 
subscript T is the counterpart in the two-channel model 
of the bound state in the single-channel model. 

The character of the other |A| — > oo bound state that 
we have denoted by the subscript A is equally obvious. 
The energy tends to oj\ — > A, the energy of the bare 
molecule, and the molecular fraction behaves as /a — > 1. 
This bound state simply represents a bare molecule that 
has decoupled from the atoms. 

The bound states at resonance A ~ make a more 
interesting talc. First, the structure of the state space 
and of the coupling between bare atoms and molecules, 
the kinematics of the problem, forces the existence of 
two bound states for the dressed molecule. Second, from 
Eq. (19) we see that the fraction of bare molecules f± in 
both bound states tends to zero when the atom-molecule 
coupling oc 5 vanishes, and to 1/2 when the atom- 
molecule coupling is strong. In dimensional quantities 
the borderline between the two cases is approximately at 
£ ~ Qp, when the strength of the atom- molecule cou- 
pling is comparable to the width of the continuum band 
of the dissociated states of the molecule. 

In the limit of strong atom-molecule coupling, |3| — > 
oo, the width of the continuum band is negligible. For- 
mally, cosg = in Eq. (8) and an effective two-level 
system for the amplitude /3 and the collective amplitude 
^2 q A q emerges. The association resonance then splits 
these two states apart. All told, the corresponding bound 
states represent association in a system that behaves as 
if there were no kinetic energy for the atoms. On the 
other hand, for weak atom-molecule coupling, 3 — > 0, 
the bound states are already far-detuned from the associ- 
ation resonance as a result of the width of the continuum 
band, and are effectively described by the single-channel 
theory. The two coexisting bound states we have de- 
noted by ± then mean that the usual bound state for 
an attractive atom-atom interaction and the repulsively 
bound state are present simultaneously. 

The two bound states also provide unexpected insights 
into modeling of Feshbach resonances. In molecular 
physics it is customary to think of Feshbach resonances 
in terms of multiple channels, whereas single-channel pic- 
tures are the norm in condensed-matter physics. Oddly 
enough, there seems to be little difference [19] between 
the predictions from single- and two-channel theories in 
common experimental situations with quantum degener- 
ate gases. In contrast, one bound state in the single- 
channel description and two bound states in the two- 
channel description is a stark qualitative difference. 

The characteristic frequency scale of the lattice 
physics [1] is the recoil frequency set by the atomic mass 
m and the lattice spacing d as en = TT 2 h/2md 2 , roughly 
10 x 27rkHz; the tunneling amplitude J is typically a 
fraction thereof. The two bound states are similar in 
character and therefore presumably easiest to detect si- 
multaneously on resonance, A ~ 0. To have the bound 



states well separated from the dissociation continuum 
we would like to have £ > Op ~ J. On the other 
hand, the resonance will overwhelm the lattice physics 
if £ S> e R , and the Bosc- Hubbard model itself may need 
to be amended [20, 21]. Overall, it appears that the 
best experimental parameters are in the neighborhood of 
£ ~ cr. Our question is, are atom-molecule couplings of 
this order available in practice? 

We proceed along the lines of Refs. [22] and [23]. Wc 
write the atom-molecule coupling as £ = aj 1 / 4 ^ 3 / 4 , where 
uj characterizes the free-space molecular physics and f2 
the lattice physics. Assuming that the Wannicr functions 
for atoms and molecules at each lattice site are ground 
states of a three-dimensional harmonic oscillator with fre- 
quencies u>i, we have = (lji^uj?,) 1 ^ /2tt. The charac- 
teristic order of magnitude of the lattice contribution is 
then of the order of the recoil frequency and, in turn, so 
should be the free-space contribution, which calls for a 
weak association resonance. 

For the Feshbach resonance, the free-space coupling for 
a given resonance is fixed once and for all. The remain- 
ing molecular frequency is uj — Att^ttlcl^ A^ A 2 fh? , where 
the background scattering length is ab and the magnetic 
field width of the resonance is Ab, while the difference 
between the magnetic moments of a bare molecule and 
two bare atoms is A^ . To achieve uj ~ cr requires a nar- 
row Feshbach resonance. Potential candidates include 
the 853 G resonance in Na [24] (Q ~ 1.6 x 27rkHz), the 
911 G resonance in 87 Rb [25] (31 x 2ttHz), and the 20 G 
resonance in 133 Cs [26] (75 x 27rHz). 

In photoassociation the atom-molecule coupling £ is 
adjustable according to the laser intensity. Borrow- 
ing from Ref. [23], we have uj = Ait 2 mL\pY 2 /H, where 
Lpa = rthKJ (47r) is determined from the low-intensity 
rate constant K oc /, and the natural linewidth of 
the molecular state is T. The broad natural linewidths 
complicate photoassociation in alkali metals [6], but the 
alkali-earth metals have narrow linewidths and, in partic- 
ular, the 192 MHz line in 174 Yb [7] is a ready candidate 
for a laser intensity of the order of 1 mW/cm 2 , so that 
Q ~ 63 x 2ttHz. 

We have discussed the stationary states of two bosons 
in a one-dimensional optical lattice within the two- 
channel Bose-Hubbard model. The particular result that 
there are two bound states, one of them effectively a 
dimer bound by attractive interactions and the other by 
repulsive interactions, provides a dramatic example of 
molecular physics in a lattice with no counterpart in free 
space. Wc have identified several Feshbach and photoas- 
sociation systems for which this unusual situation might 
be observable. 
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